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SOLUTION  AND  SENSITIVITY  ANALYSIS 
OE  RESDYN  MODEL  la  USING  SENSUMT 

by 

Anthony  V.  Fiacco 
Frank  A.  Amodeo 


1 .  Introduction 

Presented  in  this  paper  are  the  results  obtained  by  using 
SENSUMT  to  solve  "Model  la,"  a  resource  allocation  model  described 
in  [1]  and  [2].  Additional  background  can  be  obtained  from  [3]  and 
[4],  and  a  general  classification  scheme  and  model  designation  are 
given  in  [5].  Model  la  denotes  the  model  that  considers  only  maintenance 
allocation.  This  work  differs  from  that  reported  in  [2],  in  that  SENSUMT 
is  used  instead  of  LINDO,  and  the  original  objective  function  is  used 
rather  than  a  piecewise  linear  approximation.  The  use  of  SENSUMT  also 
allows  sensitivity  results  to  be  obtained.  An  extension  of  Model  la, 
incorporating  manpower  allocation,  is  reported  in  [6]  but  not  pursued 
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2.  Statement  of  the  Probl 


em 


Using  the  notation  and  simplified  equations  from  [1]  and  the  def¬ 
inition  of  the  objective  function  given  in  [2],  the  problem  can  be 
formulated  as  a  standard  nonlinear  program: 


maximize  f(x,a,  m,a0)  =  (1,2, 1,1,2) 
(x,a,m) 


+  aQm  )  (1  -  m1/a1) 
(a2  +  ct0m2)  (1  -  m2/a2) 

(a3  +  aom3^  O  _  m3^ a3^ 

(a4  +  aom4}  (1  "  m4/a4) 

(a5  +  a0m5)  "  ^5/a5) 


ao  ^ 


a„ 


a.r 


0 ir 


a„ 


subject  to: 

(1)  _a3a1-_1  +  at  +  xt  =  B0(-1  +  ’ 

“2 

(2)  -o^at_i  “  [cl y  ~~  a^CL-j ]  m^^  +  m^  4*  x^_  =  0 

(3)  -a4at^1  -  [1  -  a4ay]  m^  +  xfc  <  0  , 


t  =  1,2, . . . ,5 
t  =  1, 2 , . . •  ,5 
t  =  1,2,. ..,5 


The  variables  x^_  ,  m^_  and  a^_  denote  the  maintenance  supplied, 

the  maintenance  backlog,  and  the  asset  value  at  time  t  ,  respectively. 
The  model  parameters  are  as  follows: 

=  0.03  (yearly  proportional  budget  growth) 

&2  =  2  (ratio  of  maintenance  backlog  on  retiring  ships  to 
maintenance  average) 

a3  =  29/30  (fraction  of  assets  remaining  each  year) 

&4  =  0.04  (fraction  of  assets  requiring  annual  maintenance) 
ay  =  0.50  (fraction  of  backlog  not  requiring  additional 
maintenance) 

Gig  =  58  (Beta-function  parameter  used  in  the  objective 
function)  . 
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The  initial  conditions  are  as  follows: 

Bg  =  8  (budget) 

ag  =  100  (asset  value) 

mg  =  1  (maintenance) 

3 .  Elimination  of  the  Inequality  Constraints 

Subtracting  Equation  (2)  from  Inequality  (3)  yields 
a2 

[a3  -  1]  mt_x  -  mt  <  0  ,  t  =  1,2,... ,5  . 

From  the  definition  of  these  parameters,  a„  <_  1  and  >_  0  .  There¬ 
of  3  2  ” 

fore,  [oig  -  1]  <_  0  .  Since  m^_  >_  0  and  m^__^  5  Inequality  (3) 

is  redundant  and  can  be  eliminated. 

4 <  Attempt  at  an  Analytic  Calculation  of  Solution 

of  RESDYN  Model  la 

For  a  time  frame  of  5  years,  the  problem  has  15  variables 

^mt5  at5  an<^  xt5  ^  =  1*2,..., 5)  and  10  equality  constraints  (plus 
the  nonnegativity  constraints).  It  is  theoretically  possible  to  use 

the  equalities  to  eliminate  10  variables,  and  to  formulate  the  problem 
in  terms  of  5  variables,  subject  to  the  15  nonnegativity  constraints. 
This  procedure  was  attempted  but  proved  analytically  intractable,  due 
to  the  complicated  expressions  that  resulted.  Therefore,  another  ap¬ 
proach  was  attempted. 

From  the  results  of  [2],  it  is  concluded  that  m^  =  0, 
t  -  1,2,..., 5  in  the  given  solution.  Assuming  this,  we  may  solve  for 
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xt  and  at>  t=l,2,...,5,  and  check  whether  the  Karush-Kuhn- Tucker 
optimality  conditions  hold  at  the  resulting  point  and  hence  whether 

this  point  is  a  candidate  solution  of  the  problem  as  formulated  above, 
(These  conditions  are  necessary  at  a  solution  point  since  the  constraints 
are  linear.)  If  these  conditions  hold,  then  we  may  conclude  that  the 
problem  may  have  been  solved  by  the  given  point,  for  the  given  set  of 
parameter  values.  Solving  Equations  (1)  and  (2),  under  the  assumption 

that  m^_  =  0,  t  =  1,2,  ...5  ,  yields  the  following.  For  notational 

a2 

convenience,  let  [a^  -  a^a_J  “  &  • 


For  t  =  1,...,5  ,  Equations  (1)  and  (2)  imply: 


X1 

— 

k  m0  +  a4aQ 

— 

4.91444 

al 

= 

B„(l  +  a1)  +  a,a0  - 

X1  = 

99.99223 

X2 

= 

Vl 

= 

3.99969 

a2 

= 

2 

Bq(1  +  a1)  +  a3a1  - 

-  x2  = 

101.14667 

X3 

= 

a4a2 

= 

4.04587 

a3 

= 

Bq(  1  +  a-j^)3  +  a3a2 

-  x3  = 

102.47106 

x4 

= 

a4a3 

= 

4.09884 

a4 

= 

4 

Bg (1  +  dj)  +  a3a3  ■ 

-  x4  - 

103.96059 

x5 

= 

a4a4 

= 

4.15842 

a5 

= 

Bq(1  +  a1)5  +  a3a4  - 

-  x5  = 

105.61101 

Substituting  our  calculated  candidate  solution,  the  value  of  the  objective 
function  is  f (x,a,m,a0)  =  (1 , 2,1,1, 2) ^ ,a2,a3,a4 ,a5) T  =  719.939. 
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It  is  interesting  to  compare  our  calculated  solution  to  the 
solution  given  in  [2],  both  of  which  are  summarized  in  the  last  two 
columns  of  Table  1.  With  the  exception  of  a^  and  the  value  of  the 
objective  function,  the  results  are  very  close.  The  objective  function 
value  of  [2]  is  higher.  We  find  that  the  discrepancy  in  the  objective 
function  is  due  mainly  to  the  different  values  obtained  for  a^_  . 
Substituting  the  variable  values  of  [2]  (x^  =  4.158567,  a^  =  103.964188, 

and  a^  =  109.772942)  into  Equation  (1)  for  t  =  5  ,  we  find  that 
the  solution  of  [2]  violates  this  equality  constraint.  Hence, 
as  given  in  [2]  is  apparently  incorrect,  and  the  solution  point 
obtained  in  [2]  is  not  feasible. 

To  check  our  calculated  solution  candidate,  the  Lagrangian  of 
the  problem  is  formulated.  A  change  in  notation  is  introduced  at  this 
point,  for  simplicity:  a^  =  Xg,  a2  =  Xy,...,a^  =  x1Q,  =  x^, 

m2  ^  ~  X15  •  Also,  since  we  are  maximizing  the  objective 

function,  the  negative  of  the  objective  function  is  used  in  the 
Lagrangian  and  in  SENSUMT,  and  the  problem  is  treated  as  one  of  mini¬ 
mization.  Since,  by  assumption,  only  m  =  0,  t=l,2,...,5,  only 

the  nonnegativity  constraints  associated  with  these  variables  are 
included . 
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TABLE  1 

PROBLEM  INITIAL  CANDIDATE  SOLUTIONS 


Hand-Calculated 
Numerical  Solu¬ 
tion 

SENSUMT 

§TM-69340/83 
(Reference  [2]) 

Calculated 
with  Assump¬ 
tion  m  =  0 

f+ 

720.050 

720.050 

728.280 

719.939 

X1 

4.88272 

4.88115 

4.914398 

4.91444 

x2 

3.99920 

3.99936 

3.999825 

3.99969 

x3 

4.04369 

4.04303 

4.046005 

4.04587 

x4 

4.09731 

4.09716 

4.098984 

4.09884 

X5 

4.15730 

4.15695 

4.158569 

4.15842 

al 

= 

x6 

100.02395 

100.025 

99.996 

99.992 

a2 

= 

X1 

101.17781 

101.179 

101.150 

101.147 

a3 

= 

x8 

102.50335 

102.505 

102.475 

102.471 

a4 

= 

x9 

103.99334 

103.995 

103.964 

103.961 

a5 

= 

x10 

105.64379 

105.646 

109.773+ 

105.611 

ml 

= 

i — 1 
i — 1 

X 

.03172 

.0332867 

0 

0 

m2 

= 

X12 

.03076 

.0320911 

0 

0 

m3 

x13 

.03155 

.0334766 

0 

0 

m4 

= 

x14 

.03168 

.0336642 

0 

0 

m5 

= 

x15 

.03140 

.0336430 

0 

0 

W1 

5.98348 

4.42901 

W2 

5.37768 

7.53408 

W3 

3 . 64464 

-3.56954 

W4 

2.85376 

.044104 

W5 

2.00030 

5.92848 

W6 

-5.98348 

-9.07268 

W7 

-5.37768 

-6.39773 

W8 

-3.64464 

-.802994 

W9 

-2.85376 

8.71874 

W10 

-2.00030 

-6.74342 

tSee  Section  4  above  concerning  the  value  of  a,-  . 

§Point  obtained  using  a  piecewise  linear  approximation  of  objective 
function.  With  a,.  corrected  (see  Section  4) ,  this  point  is  probably 
optimal  for  this  approximation.  It  is  not  optimal  for  the  original  problem. 
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The  resulting  Lagrangian  is 


L  (x,u,w,e)  =  (-1,-2, -1,-1, -2) 


(xfi  +a0xu)(l  -  x^/xj 


ar 


11  6 


(xy  +  a0x12)(l  -  x12/xy) 

(x8  +  a0x13)(l  -  x13/x8) 

(x9  +  a0x14)(1  “  x14/x9} 


ar 


ar 


(x10  +  a0x15)ci  '  x15/x10} 


a„ 


U1XH  U2X12  ”  U3X13  "  U4X14  "  U5X15 


+  w1  [x3 


+  X6  -  B0(l  +  op  -  a3a0] 


+  «2  [x2  -  CIjXj  +  X?  -  Bq(1  +  o3)  ] 

+w3  [x3-a3x7+x8  bq(1  +  c^) 2] 

+  w4  [x4-a3x8+x9  “  Bq(1  +  a1)3] 
+  w5  [x5  -  a3x9  +  x1Q  -  B0(l  +  a^5] 


+  w*  [x. 


+  X. 


11 


a4a0  “  km0] 


+  w7  [x2  -  a4x6  -  kxn  +  x12] 


+  WR  [x,  -  a4x7  -  lac-,0  + 


12  13 


+  wg  [x4  -  a4x8  -  kx13  +  x14l 


+  w10  [x5  a4x9  kx14  +  x15] 

where  £  =  (o^,  a2,  a3,  a^,  a?,  BQ,  m0>  aQ,  aQ) 


Taking  the  gradient  of  the  Lagrangian  with  respect  to  the  15  variables  gives 
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0 


0 


(-1) 

(-2) 

(-1) 

(-1) 

(-2) 


2  2 

1  +  (otg-l)  (x1;l/x6)  +  aQ  (xi;L/x6) 


(l-(xll/x6)) 


a0-1 


1  +  (o£q-1)  (x15/x10)  +  a0  ^X15^X10^ 


(l-(x15/xl0) 


ao-1 


(+1)  ao^ao+1^  ^xl]/x6^  (1-(xii/x6) 


L 


(+2) 

(+1) 

C+l) 


a  -i 

(+2)  aQ(a0+l)  (x15/x10^  (l-(x15/x10^ 


0 

W1 

+ 

W6 

• 

W2 

+ 

W7 

• 

W3 

+ 

W8 

• 

W4 

+ 

W9 

0 

W5 

+ 

W10 

0 

W1 

- 

a3w2  ‘  a4W7 

• 

W2 

— 

a3w3  -  a4wg 

• 

+ 

W3 

— 

a3w4  "  a4w9 

• 

w4 

— 

a3w5  “  a4w10 

0 

w5 

U1 

W6 

■  k 

W7 

U2 

W7 

-  k 

W8 

U3 

W8 

■  k 

w9 

u4 

w9 

■  k 

W10 

U5 

1 

M 

O 
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Let  xn  =  x12  =  x13  =  x14  =  x15  =  0  .  Then 


The  gradient  of  the  Lagrangian  is  set  equal  to  zero  to  see  if 
the  Karush-Kuhn- Tucker  first  order  necessary  conditions  hold.  From 
row  10,  w5  =  2-.  From  row  5,  w1Q  =  -2  .  From  row  15,  u  =  -2  . 
Since  the  multipliers  associated  with  inequality  constraints  must  be 
nonnegative  for  the  point  to  be  a  minimizer  (maximizer  of  the  orig¬ 
inal  problem) ,  the  calculated  candidate  solution  point  is  not  even  a 
local  minimum.  Hence,  the  variables  cannot  be  equal  to 

0  at  a  solution  of  RESDYN  Model  la. 

It  should  be  noted  in  passing  that  [2]  was  obtained  by  using 
a  piecewise  linear  approximation  of  the  objective  function.  If  the 
gradient  of  the  Lagrangian  is  reformulated  using  this  approximation, 
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it  may  be  shown  that  the  first  order  necessary  conditions  are  satis¬ 
fied  at  our  calculated  solution  point. 

Returning  to  our  analysis  of  the  original  problem,  since  none 
of  the  nonnegativity  constraints  are  binding,  it  follows  that  all  the 
u^  =  0  .  We  now  find  that  the  gradient  of  the  Lagrangian  can  be  set 
equal  to  zero  and  that  we  can  solve  analytically  for  the  multipliers 
Wj  and  the  ratios  •  For  notational  con¬ 
venience,  define  z^  =  x1]_/x5»  z2  =  ’  z5  =  x15//x10  * 


From  row  5: 
From  row  10: 

From  row  15: 

Therefore, 

so  z  ^  = 


W5  “W10 


(-2) 


2  2 

1  +  (aQ-l)z5  +  aoz5 


(l-z5) 


V1 


+  w<_  =  0 


V1 

(+2)  a0(l+aQ) z5(l-z5)  +  =  0 


1  +  (a0-l)z5  + 

2  2 


a0  1  an  1 

(l-z5)  «  a0(l+a0)z5(l-z5) 


’5=1  or  1  +  (aQ-l)z5  +  aQz5  -  a0(l+a0)z5  =  0 
yielding 

Using  the  formula  for  the  solution  of  a  quadratic  equation  yields 


2  2  2 

aQz5  -  (aQ+l)z5  +1=0. 


z5  = 


(<Xq+1)  ±  (a^-l) 


2ar 


and  hence  z,.  =  1  and  =  w^q  =  0  or  z,.  =  1/cIq  .  It  may  be  argued 

that  the  smaller  value  of  z±  will  reduce  f  ,  hence  =  1  is 

discarded.  Substituting  back  into  either  row  10  or  15  with 
2 

z^  =  1/cXq  yields 

W5  =  "w10  =  2(1+1/ao^1-1/ao^ 

Since  w^  =  -w^  from  row  4,  the  values  for  w<.  and  w^q  can  be  sub¬ 
stituted  into  rows  9  and  14  to  solve  for  z^  .  However,  these  equations 


2.  a0~X 
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cannot  be  solved  in  closed  form,  and  in  particular  one  cannot  as  before 

v1 

divide  by  (1-z^)  to  obtain  a  quadratic  equation.  The  equations 

were  solved  numerically,  using  Newtonrs  method,  substituting  the  nominal 
values  for  the  parameters.  Substituting  the  calculated  values  of  the 
multipliers  into  the  other  rows  (8  and  13,  then  7  and  12,  etc.),  all  the 
ratios  and  multipliers  were  solved  by  this  numerical  method.  Using  the 
numerical  values  for  z^,  and  fixing  the  parameters  at  thier 

nominal  values,  the  recursive  equations  (1)  and  (2)  of  Section  2  were 
solved  to  obtain  a  candidate  solution  point.  The  optimal  value,  can- 
didate  solution  .point.  The  optimal  value,  candidate  solution  point 
and  multipliers  resulting  from  this  method  are  given  in  Table  1. 

5 .  Solution  using  SENSUMT 

The  problem  was  solved  using  SENSUMT  running  on  a  DEC  VAX  11/780. 
The  method  selected  for  minimizing  the  unconstrained  penalty  function 
was  G.  P.  McCormickTs  modification  of  the  Fletcher-Powell  method  (i.e., 
NEX0P2  =  4)  .  The  results  are  summarized  in  Table  1.  For  comparison, 
and  although  they  are  not  feasible,  the  results  of  [2]  and  the  results 
that  follow  from  the  closed  form  solution  of  the  equations  with  the 
assumption  that  mt,==  t  =  1,2,, ..,5  ,  are  also  presented,  as  noted 
earlier. 

The  solution  point  and  estimates  of  the  multipliers  were  compared 
to  the  results  of  the  method  described  above  in  Section  4  (forming  the 
gradient  of  the  Lagrangian  with  respect  to  the  15  variables  and  set¬ 
ting  this  to  zero) * 


11 


T-485 


The  SENSUMT  estimates  of  the  multipliers  do  not  appear  to  be 

accurate.  These  estimates  are  obtained  from  the  ratio  w.  =  2h./p  , 

3  3 

where  p  >  0  is  the  penalty  function  algorithm  parameter  set  by  the 
user.  Since  both  the  constraint  value  h.  and  p  go  to  zero  as  the 
program  iterates,  the  resulting  ratio  of  two  small  numbers  may  not  be 
accurate.  The  accuracy  of  the  solution  point  is  much  better  than  that 
of  the  multipliers.  All  the  variables  agree  to  the  third  place  after 
the  decimal  point.  The  optimal  value  (720.050)  is  in  agreement  to  the 
accuracy  kept  in  the  numerical  solution. 

6 .  Sensitivity  Analysis 

The  use  of  SENSUMT  allows  sensitivity  results  to  be  obtained. 

The  sensitivity  of  the  optimal  value,  solution  point  and  multipliers 

m 

were  obtained  for  nine  parameters:  e  =  (o^,  a^,  a^,  BQ,  aQ ,  mQ,  aQ) 

where  aQ  is  the  parameter  in  the  objective  function.  The  gradients 
of  the  optimal  value  function  f  ,  the  solution  point,  and  Lagrange 
multipliers  are  given  in  Tables  2,  3  and  4,  respectively.  Note  that 
for  the  sake  of  the  sensitivity  analysis,  the  initial  conditions 
(Bq j  3q  and  iHq)  are  treated  as  parameters. 

The  above  gradients  represent  a  rate  of  change  in  the  variable 
with  respect  to  each  parameter.  These  gradient  values  are  valid  only 
at  the  nominal  parameter  values,  and  can  change  significantly  as  a 
parameter  is  varied  from  its  nominal  value.  Also,  some  parameters  can 
vary  over  a  much  wider  range  than  others.  For  example,  a  parameter 
such  as  ap  (initial  asset  value),  with  a  nominal  value  of  100,  can 
change  over  a  larger  range  than  (proportional  budget  growth) ,  which 

has  a  nominal  value  of  0.03  . 
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TABLE  2 

SENSITIVITY  OF  f  (e) 


419,141 

ai 

-10,1688 

a2 

2052.00 

, 

a3 

-1836.59 

aA 

where  £  = 

-12.8398 

a7 

21.3820 

Bo 

5.54479 

ao 

-5.47166 

mo 

26.5243 

JV 
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TABLE  3 

SENSITIVITY  OF  SOLUTION  POINT 


VTx  = 
s 


vta  = 


Vjm  = 


2. 2E-5 

,32QQ 

,9558 

1,904 

3,163 

-235,5 

197,3 

350,8 

-563,4 

242.2 

5597, 

120.5 

-25145. 

27484, 

-8041. 

3795. 

-3030. 

-5316. 

8891. 

-3733. 

-297.3 

249.1 

443.0 

-711.4 

305.8 

3.0E-6 

•  4.1E-2 

8.1E-2 

.1184 

.1547 

4. OE-2 

3.7E-2 

3.4E-2 

3.2E-2 

2.9E-2 

.9144 

-3.7E-2 

-3.4E-2 

-3.1E-2 

-2.9E-2 

43.95 

36.02 

-30.78 

2.044 

44.18 

8.000 

23.89 

47.60 

79.08 

118.3 

235.5 

30.37 

-321.5 

252.6 

2.048 

-5497. 

-5334. 

20090. 

-7962. 

448.0 

-3795. 

-638.2 

4699. 

-4349. 

-471.2 

297.3 

38.35 

-405.9 

319.0 

2.586 

1.030 

2.015 

2.960 

3.869 

4.744 

.9267 

.8587 

.7957 

.7374 

.6833 

-.9144 

-.8474 

-.7852 

-.7277 

-.6743 

_  -43.95 

78550 

-45.10 

-45.64 

-88.30 

-2.2E-5 

-3.9E-5 

-5.9E-5 

-7.2E-5 

2.0E-4_ 

235.5 

27.46 

-324.5 

253.8 

1.6E-2 

-5595. 

-5456. 

19942. 

-8445. 

-.5265 

-3695. 

-40Q.8 

5025. 

-4006. 

-.2443 

297.3 

34.67 

-409.7 

320.5 

2. OE-2 

-3.0E-6 

-5.4E-6 

-2.4E-6 

-5.1E-6 

9.7E-6 

-7.1E-7 

2.1E-7 

2.2E-7 

2.8E-7 

1.9E-6 

4.5E-5 

-2.0E-7 

-2.1E-7 

-2.8E-7 

-1.8E-6 

-43.95 

-77.96 

-43.65 

-43.76 

-86.02 

Recall  that 


3..^  x6 ,  •  •  •  ,  a^  —  x 


10  ’  ml  xn » 


,m5  x15 
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TABLE  4 

SENSITIVITY  OF  LAGRANGE  MULTIPLIERS 


-2.7E-8 

-2.8E-8 

-2.6E-8 

-1.9E-8 

-1.1E-8 

-5.4E-7 

-6.0E-8 

8 . 2E-7 

-5.4E-7 

-1.9E-9 

8.982 

1.562 

5.538 

5.928 

1.2E-7 

-5.021 

1.486 

2.470 

-6.743 

1.6E-7 

-7.0E-7 

-8.2E-8 

1.0E-6 

-6.6E-7 

3.4E-9 

-6.6E-9 

-6.8E-9 

-3.8E-9 

-3.6E-9 

-2.9E-10 

-5.8E-10 

-4.4E-10 

-2.9E-10 

-1.7E-10 

-5.8E-11 

5.8E-10 

4 . 4E-10 

2.9E-10 

1.5E-10 

2.9E-11 

3.5E-5 

3.1E-5 

2.1E-5 

1.6E-5 

1.1E-5 

2.8E-8 

9 . 4E-9 

-3.0E-8 

-8.7E-8 

-1.7E-7 

1.0E-5 

1.2E-5 

-2.1E-5 

3.2E-5 

-1.3E-5 

-8.982 

-1.562 

-5.536 

-5.930 

4.4E-4 

5.021 

-1.486 

-2.470 

6.743 

2.1E-4 

1.3E-5 

-1.5E-5 

-2.7E-5 

4.1E-5 

-1.7E-5 

6 . 6E-9 

4.4E-9 

-7.0E-10 

-3.0E-9 

-8.2E-9 

-9.9E-10 

-1.7E-9 

-1.7E-9 

-1.6E-9 

-1.6E-9 

-3.7E-8 

1.7E-9 

1.7E-9 

1.6E-9 

1.6E-9 

-3.7E-5 

-3.3E-5 

-1.9E-5 

-1.7E-5 

-1.4E-5 

1  T 

W  —  y  .  ,  .  ?W^) 

W2  =  (w6,...,wi0)T 
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In  order  to  identify  those  parameters  that  have  the  greatest 
impact  on  the  solution  of  the  problem,  it  is  useful  to  calculate  for 
each  parameter  the  change  in  the  optimal  value  resulting  from  a  change 
of  a  fixed  proportion  of  the  range  associated  with  the  parameter.  An 
estimate  of  the  range  for  each  parameter  was  obtained  from  [7].  These 
ranges  are  given  in  Table  5.  For  each  parameter,  an  estimate  of  the 
change  in  the  optimal  value  due  to  an  increase  in  the  parameter  of  1% 
of  the  range  (while  all  other  parameters  are  held  at  their  nominal 
values)  is  calculated.  This  calculation  is  based  on  an  extrapolation 
of  the  first  order  sensitivity  estimates  at  the  solution,  i.e., 

•k 

v£f  *  Since  this  extrapolation  is  linear,  the  change  in  the  optimal 
value  function  is  proportional  to  the  proportion  of  the  range  selected. 
Although  the  proportion  chosen  is  1%,  the  extrapolation  could  be  per- 
formed  for  any  proportion.  These  changes  in  f  due  to  each  parameter 
are  given  in  Table  6. 

It  is  useful  to  rank  the  parameters  in  order  of  their  impact  on 
the  optimal  value,  using  the  absolute  value  of  the  elements  of  the 
gradient  and  the  changes  of  the  optimal  value  due  to  a  change  of  a 
fixed  proportion  of  the  range  of  each  parameter.  For  the  gradient,  the 
order  of  sensitivity  to  the  parameters  (from  greatest  to  least  sensi¬ 
tivity)  is  a^,  a4,  a^,  a^,  Bq,  a^,  ag,  and  m^  .  Notice  the 

dominance  of  ex ^  and  on  the  optimal  value,  and  also  on  the  solu¬ 

tion  point  and  multipliers  (from  Tables  3  and  4).  To  a  lesser  extent. 
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TABLE  5 

RANGES  ASSOCIATED  WITH  THE  PARAMETERS 


Parameter 

Nominal 

Value 

Lower 

Bound 

Upper 

Bound 

Rang' 

ai 

0.03 

-0.05 

0.15 

0.20 

a-2 

2.0 

1.0 

4.0 

3.0 

a3 

0.967 

0.75 

1.0 

0.25 

a4 

0.04 

0.01 

0.10 

0.09 

a7 

0.50 

0 

1.0 

1.0 

Bo 

8 

5 

20 

15 

ao 

100 

80 

300 

220 

mo 

1 

1 

10 

9 

ao 

58 

5 

60 

55 

17 


T-485 


TABLE  6 


CHANGES 

IN  f*(e)  FOR  AN 

INCREASE  OF 

1% 

OF  THE  RANGE  FOR  EACH 

PARAMETER 

Parameter 

V*(£) 

* 

Af 

*  * 
f  +  Af 

ai 

419.141 

.8383 

720.8883 

a2 

-10.1688 

-.3051 

719.7449 

a3 

2052.00 

5.1300 

725.1800 

a4 

-1836.59 

-1.6529 

718.3971 

a7 

-12.8398 

-.1284 

719.9216 

R0 

21.3820 

3.2073 

723.2573 

a0 

5.54479 

12.1985 

732.2485 

mo 

-5.47166 

-.4924 

719.5576 

ao 

26.5243 

14.5884 

734.6384 
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the  optimal  value  is  sensitive  to  ,  while  the  solution  point  is 

sensitive  to  and  a ^  .  All  other  sensitivities  are  quite  low. 

for  the  changes  in  the  optimal  value  due  to  a  range  proportional  change 
in  each  parameter,  the  order  of  sensitivity  is  aQ ,  aQ ,  a3 ,  BQ ,  , 

iHq,  and  oty  .  Notice  how  the  ranking  changes  due  to  the  variance 

in  the  ranges.  Now  and  a^  are  the  most  critical  parameters  in 

the  optimal  value. 

In  order  to  get  an  estimate  of  lower  and  upper  bounds  on  the 
optimal  value,  the  SENSUMT  program  was  run  with  the  parameters  set  to 
their  extreme  values.  The  particular  extreme  value  (upper  or  lower 
bound)  depended  on  the  sign  of  the  corresponding  element  in  the  gradient. 

X 

The  estimate  for  a  lower  bound  on  f  (e)  is  248.844.  However  at 
this  solution  point  the  sign  of  the  gradient  changes  for  a  9  c^,  and 
QL-j  compared  to  the  gradient  at  the  nominal  solution  point.  Therefore 
it  is  possible  to  obtain  an  even  lower  optimal  value  at  some  inter¬ 
mediate  values  of  these  parameters.  The  estimate  of  the  upper  bound 
is  2647.162.  All  elements  of  the  gradient  have  the  same  sign  as  the 
gradient  at  the  nominal  solution.  Therefore,  it  is  plausible  that 
this  estimate  of  the  upper  bound  is  valid. 

Calculation  of  parametric  optimal  value  bounds  for  this 
nonconvex  program  and  analyses  for  extensions  of  this  model  will  be 
a  subject  of  future  research. 


-  19  - 


T-485 


REFERENCES 

[1]  FALK,  J.  E.  On  solving  toy  model  la.  Technical  Memorandum 

Serial  TM-69310/83,  16  February  1983,  Institute  for 
Management  Science  and  Engineering,  The  George  Washington 
University. 

[2]  FALK,  J.  E.  and  G.  P.  McCORMICK.  Solution  of  model  la  with  beta 

data.  Technical  Memorandum  Serial  TM-69340/83,  7  March  1983, 
Institute  for  Management  Science  and  Engineering,  The  George 
Washington  University. 

[3]  FALK,  J,  E.  and  G.  P.  McCORMICK.  Solution  of  problem  la  with 

different  beta  data.  Interdepartmental  Memorandum  Serial 
69351,  11  March  1983,  Institute  for  Management  Science  and 
Engineering,  The  George  Washington  University. 

[4]  CLARK,  R.  Resource  dynamics  project  status.  Technical  Memo¬ 

randum  Serial  TM-69357/83,  17  March  1983,  Institute  for 
Management  Science  and  Engineering,  The  George  Washington 
University. 

[5]  CLARK,  R.  RESDYN  optimization  models  and  model  designation. 

Technical  Memorandum  Serial  TM-69455/83,  23  May  1983, 
Institute  for  Management  Science  and  Engineering,  The  George 
Washington  University. 

[6]  FALK,  J.E.,  and  G.  P.  McCORMICK.  On  solving  a  dynamic  optimal 

budgeting  problem.  Technical  Paper  Serial  T-474/83,  9  June 
1983,  Institute  for  Management  Science  and  Engineering, 

The  George  Washington  University. 

[7]  CLARK,  R.  Informal  correspondence,  1984. 


20 


